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We present a microscopic quantum analysis for rotational constants of the OCS-H2 complex in 
helium droplets using the local two-fluid theory in conjunction with path integral Monte Carlo 
simulations. Rotational constants are derived from effective moments of inertia calculated assuming 
that motion of the H2 molecule and the local non-superfluid helium density is rigidly coupled to the 
molecular rotation of OCS and employing path integral methods to sample the corresponding H2 
and helium densities. The rigid coupling assumption for H2-OCS is calibrated by comparison with 
exact calculations of the free OCS-H2 complex. The presence of the H2 molecule is found to induce a 
small local non-superfluid helium density in the second solvation shell which makes a non-negligible 
contribution to the moment of inertia of the complex in helium. The resulting moments of inertia for 
the OCS-H2 complex embedded in a cluster of 63 helium atoms are found to be in good agreement 
with experimentally measured values in large helium droplets. Implications for analysis of rotational 
constants of larger complexes of OCS with multiple H2 molecules in helium are discussed. 

PACS numbers: 36.40.-c, 36.40.Mr, 67.40.-w, 67.40.Yv 
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I. INTRODUCTION 

Helium droplets are known to provide an attractively 
convenient environment for synthesis of weakly bound, 
metastable, or other unusual complexesi. The low tem- 
perature of these droplets (~ 0.15 — 0.4 K), their liquid 
nature, and their weak interaction with impurity species 
render them ideal matrix hosts for weakly bound com- 
plexes. Formation of such complexes is facilitated by the 
now standard pick-up technique 2 in which foreign species 
are absorbed by the droplet in a pick-up chamber. Se- 
quential use of such pick-up allows assembly of a variety 
of complexes and aggregates of atoms and molecules. In 
recent years this synthetic potential of helium clusters 
has been demonstrated with synthesis of water clusters^ 
of metal clusters? 4 *^, of metal-adsorbate clusters^ and of 
chains of HCN, 7 all in the ultra-cold environment of a 
helium droplet. In several of these instances high reso- 
lution spectroscopy of the embedded clusters has shown 
that metastable isomers of conformations not seen for 
the corresponding clusters in the gas phase are formed 
in helium. This indicates that although weak, the effect 
of helium solvation is non-negligible in determining the 
structure of some complexes, in particular those formed 
from polar constituents. 

Complexes of molecules with molecular hydrogen con- 
stitute a special type of complex whose formation and 
behavior in helium droplets is currently of great interest. 
The possibility of finding a superfluid state of molecular 
hydrogen provides an additional motivation beyond that 
of merely forming the complexes and analyzing the effect 
of helium solvation on their structure. Molecular hydro- 
gen superfluidity is most likely to be found for the j = 
state of H2, i.e., para-H2, the lightest bosonic molecule 
of the hydrogen isotopes. We shall restrict ourselves to 



para-H2 in this work, and for simplicity refer to it as 
H2. Helium will in all cases be taken as the boson iso- 
tope, He. A number of significant experimental results 
have been obtained for OCS(H2)a/ complexes in recent 
years. Ref. showed that for M > 11 there exist spec- 
tral anomalies which are consistent with the existence of 
a superfluid state of the H2 molecules. Existence of this 
superfluid state of molecular hydrogen in a nanoscale sys- 
tem has recently been confirmed by path integral calcu- 
lations and shown to have an onset at a temperature of 
~ 0.3 K for hydrogen clusters around OCS£ Additional 
spectroscopic studies have been made for M = 1 — 8 in 
helium droplets and analyzed in terms of rigidly coupled 
models of the OCS(H2)a/ complexes i&lL For the small- 
est M — 1 complex of OCS-H2, the corresponding gas 
phase spectroscopy has recently been measured and an- 
alyzed in Ref. [l2- This smallest complex provides a ref- 
erence point for theoretical analysis of this and all larger 
clusters with M > 1 in helium. 

In this work we undertake a detailed theoretical anal- 
ysis of the M = 1 complex OCS-H2 in droplets of 4 He, 
using path integral methodology and the local two-fluid 
theory^ We make path integral Monte Carlo calcula- 
tions of the complex in clusters of N = 63 helium atoms. 
This size is large enough to provide a complete first sol- 
vation shell around the OCS molecule (~ N = 20) 14 
that has a robust and compact structure which is in- 
dependent of further increases in N. Both hydrogen 
molecule and helium atoms are described quantum me- 
chanically, with full permutation exchange symmetry of 
the helium component incorporated. We compute the 
local non-superfluid helium solvation density induced by 
the OCS-He and H2-He interactions. Comparison of this 
non-superfluid density around OCS in pure 4 He64 and 
around OCS-H2 in 4 He63 shows that the H2 molecule in- 
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duces a small non-superfluid density around itself in the 
second solvation shell of helium as well as in the first 
shell. 

We evaluate the effective moment of inertia of the OCS 
molecule assuming that the H2 molecule is rigidly coupled 
to the molecular rotation and employing the microscopic 
two-fluid theory of Ref. to evaluate the contribution 
from the helium solvation density. The assumption of 
rigid coupling of H2 to the molecular rotation is tested 
by an evaluation of the rotational constants of free (gas 
phase) OCS-H2 and by comparison of these with the ex- 
perimental values measured in Ref. 0. This compari- 
son, made with interaction potentials^ that have been 
previously calibrated for OCS-H 2 by exact bound state 
calculations,- shows that the free OCS-H 2 complex can 
be very accurately described with a rigid coupling model. 
This free complex is planar, with mass distribution cor- 
responding to an asymmetric top. Excellent agreement 
of the calculated moments of inertia with experimental 
values is found. The free complex calculations are fol- 
lowed by a local two-fluid analysis of the response of the 
solvating helium density to rotation of this OCS-H2 com- 
plex, to obtain estimates for the helium contribution to 
all three moments of inertia in helium droplets. The av- 
erage moment of inertia resulting from this analysis in- 
corporating both first and second shell complex-induced 
non-superfluid helium densities is in excellent agreement 
with the corresponding experimentally measured value 
in helium droplets, with an accuracy of ~ 4%. The cal- 
culated inertial defect indicates a non-planar mass dis- 
tribution of the total hclium-solvated complex, in agree- 
ment with the experimental measurements. Comparison 
of individual moments of inertia with their corresponding 
experimental values shows somewhat larger deviations of 
10 - 20%. The smaller accuracy of the individual mo- 
ments of inertia can be ascribed to limitations of the lo- 
cal non-superfluid estimators used here in distinguishing 
accurately between rotations around different axes. 

Such quantitative agreement of the moments of iner- 
tia and of the corresponding rotational constants with 
experimental measurements implies that the underlying 
model of an OCS molecule rotating with a rigidly coupled 
H2 molecule and a rigidly coupled local non-superfluid 
density of helium provides an accurate description for 
the low rotational states of the OCS-H2 complex that 
are accessed by spectroscopic measurements. Given that 
high quality interaction potentials calibrated by experi- 
mental measurements of the free OCS-He and OCS-H2 
complexes are employed here, we can take the agreement 
between these results for the hydrogen complex of OCS 
embedded in helium with experiment as providing further 
evidence for the accuracy of the local two-fluid model of 
quantum solvation and superfluid response to rotational 
motion around heavy molecules .i&ii 

The structure of the rest of the paper is as follows. 
In Section [n] we summarize the path integral method- 
ology used for these calculations with both helium and 
hydrogen treated quantum mechanically. In Section IIIII 



we describe the interaction potentials and present the 
various calibration calculations for the gas phase OCS- 
H2 complex that provide a validation of our assumption 
that at low energies the dynamics of the H2 molecule 
are rigidly coupled to the molecular rotation. Section llVl 
then presents the results of the full path integral calcula- 
tions and microscopic two-fluid analysis for the OCS-H2 
complex when embedded in helium clusters. We conclude 
in Section Ivl with a summary and brief discussion of im- 
plications for larger complexes of OCS(H2)m in helium. 



II. METHODOLOGY 

In this work the OCS molecule is treated as fixed 
and non-rotating. These restrictions are convenient be- 
cause of the greater computational cost associated with 
the quantum treatment of the additional H2 molecule. 
The consequences of these restrictions are not severe 
for the two-fluid analysis, as the following arguments 
indicate. Previous studies have shown that for heavy 
molecules the helium solvation density around an embed- 
ded molecule is not significantly affected by the molecule 
translation motionjiSiiS but that there is a significant de- 
crease in the angular modulation of this solvation den- 
sity when the molecular rotational degrees of freedom are 
incorporatedii^S Very recent path integral studies that 
incorporate all molecular translational and rotational de- 
grees of freedom have shown that the integrated local 
non-superfluid density is insensitive to this modulation 
and that the moment of inertia contributions from the 
local non-superfluid may be accurately estimated from 
calculations made for a non-rotating and non-translating 
OCSiSi Therefore we shall employ a non-rotating OCS 
molecule here and postpone investigation of the detailed 
effects of OCS rotational and translational motions to 
future work. 

For the interaction potentials between the three differ- 
ent particles involved (para-H 2 , He, and OCS), we use 
a sum of pair potentials consisting of H 2 -Hej22i He-He^ 
OCS-H 2 ^ and OCS-He^ terms. The para-hydrogen 
molecules are treated as spherical particles like the he- 
lium atoms. This is a well justified approximation at 
these low temperatures because of the large rotational 
constant and nearly isotropic ground-state electronic con- 
figuration of H2. Consequently, we average the three- 
dimensional ab initio OCS-H2 potential which has re- 
cently been computed by Higgins et al. with fourth- 
order M0ller-Plesset (MP4) perturbation theory^ over 
the orientation of the H2 molecule. This results in a 
two-dimensional representation of the OCS-H2 interac- 
tion that has a very similar topology to the previously 
reported MP4 OCS-He potential^ Contour plots of the 
two potentials are shown in Figure ^ It is evident that 
the OCS-H2 potential has its global minimum at approx- 
imately the same location as the OCS-He potential, but 
possesses a significantly deeper value. For OCS-H2, the 
global minimum is —208 K, located at r = 3.35 A and 
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FIG. 1: Contour plots of the H 2 -OCS and He-OCS poten- 
tials, shown as a function of the cylindrical coordinates (Z,R), 
where Z is the coordinate along the axis passing through the 
OCS molecule and R is the radial distance from the OCS 
molecular axis. The origin is set at the center of mass of the 
OCS and the molecule is oriented as O-C-S from — Z to +Z. 
Contours are shown in increments of 15 K for the H2-OCS 
and 6 K for the He-OCS potential, respectively. 



9 = 106°, where r is the distance from the OCS center of 
mass and 9 the polar angle measured from the sulfur side 
of the OCS molecular axis. For OCS-He, the global min- 
imum is -65.3 K and is located at r = 3.38 A, 9 = 108°. 
All calculations described here employ these two MP4 
potentials for the interaction of OCS with H2 and He. 
For the H 2 -He and He-He interactions we use the empiri- 
cal potential proposed by van den Bergh and Schouten^ 
and the semi-empirical potential of Aziz et al.^ respec- 
tively. 

We now have the following system Hamiltonian for the 



H2-OCS complex inside the 4 Hejv cluster: 

N N 

+ y^yk 2 -He(nM) + %2-ocs(n>) + ^Vnc-ocsin)- 

i=l i=l 

(1) 

Here mu 2 (tihc) is the mass of a hydrogen molecule (he- 
lium atom) and r*o is the position vector of the hydrogen 
molecule in a (space fixed, for a non-rotating molecule) 
frame centered on the OCS molecule. 

To study the OCS-H2 complex and its helium solvation 
structure we employ here the path integral Monte Carlo 
(PIMC) technique. This approach allows us to make the 
quantitative analysis of superfluidity in the helium envi- 
ronment that is required for the local two-fluid theory 
The main elements of the PIMC technique have been de- 
scribed in Ref. and its adaptation to molecule-doped 
helium clusters detailed in Ref. 0. We give here only a 
brief summary, with methodological details that are spe- 
cific to the OCS-H2 and its helium surroundings. The 
thermal density matrix is given by 

P {n,K';f}) = {n\e-P H \n'), (2) 

with (3 = 1/ksT, 1Z a 3(N + l)-dimensional vector, 
1Z = (fb, ■■■'Tn), and H equal to the Hamiltonian, 
Eq. J3>, for an H2 molecule and N helium atoms in the 
external field provided by the stationary OCS molecule. 
In order to incorporate the bosonic symmetry of the 4 He 
atoms, the density matrix is symmetrized by summing 
all permutations V amongst the N helium atoms: 

p B {n,n':p) = ^Y.p{n,rn'-p). (3) 

V 

Since the density matrix of a interacting quantum sys- 
tem is generally not known at a low temperature T, it 
is replaced in the Feynman path integral representation 
with a product of L higher-temperature density matrices, 
resulting in the discrete time path integral: 

p(1Z, VIZ'; (3) = J ■ ■ J dKidR.2 ■ ■ ■ dK L -i 

piTZ, Hi;t)p(Ki ,1Z 2 ;t) ■ ■ ■ p(lZ L -i,Vn'; t). (4) 

Here t = fi/L constitutes the imaginary time step 
dehning this discrete representation of the path integral. 
At sufficiently high temperatures LT, or equivalently, 
at sufficiently small time steps r, the density matrix 
p(1Z,1Z';t) can be approximated by a product of the 
free particle propagator and an interaction term. For 
the spherical He-He and H2-He interactions, the pair- 
product form of the exact two-body density matrices can 
be used, and a matrix squaring methodology employed 
to achieve accurate grid representation of these^ 
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The anisotropic He-OCS and H2-OCS interactions are 
treated here within the primitive approximation*^ From 
our previous study of OCS-doped He at clusters, we have 
found that a time step of r _1 = 80K is required to get 
converged helium density profiles around the molecule 
with this approximation for the high-temperature 
density matrix, 2 ^ In order to make a simultaneous calcu- 
lation of the sum over ./V-particle permutations in Eq. @ 
together with the multi-dimensional integration over the 
position coordinates in the discrete time path integral, 
Eq. @, we employ a stochastic process which sam- 
ples the discrete paths {1Z, Hi,1li, . . . ,1Zl-i,VTZ'} 
with probability density proportional to 
p(K, Tlr,T)p(Ki,n 2 ; r) • • • p(K L -i,VR,'; r). This 
sampling is performed by the generalized Metropolis 
sampling algorithm of Pollock and CeperleyiSi A detailed 
description of this algorithm is provided in Refs. 28, 25 
and [l^. The thermal average of any observable O can 
then be determined by computing an arithmetic average 
of (1Z'\0\TZ) over the paths sampled. 

One of the advantages of using the path integral ap- 
proach is that it allows a quantitative estimate to be 
made of the superfluidity of bosonic systems such as 
4 He and para-H2 clusters as a function of temperature. 
Within the Feynman path integral analysis, the global 
superfluid fraction that is defined as the linear response 
of the bosonic system to classical rotation of its bound- 
aries can be evaluated by an estimator written in terms 
of the projected area of the Feynman paths 2 ^: 

Here Ai is the area of a Feynman path projected onto a 
plane perpendicular to the axis Xi, and I cl is the clas- 
sical moment of inertia tensor. This estimator yields a 
non-negligible value only when exchange-coupled paths 
are comparable to the size of the system. One problem 
that arises when applying this global superfluid estimator 
to analysis of doped helium clusters is that it does not 
give information about the local perturbation of super- 
fluidity due to the presence of an impurity. However, it 
is precisely the local perturbation of the superfluid that 
determines the dynamic response of the helium to the 
rotation of an embedded moleculeji^iii so that a local 
analysis of superfluidity in the solvating helium density 
is required. Another possible issue with this estimator is 
that it provides a measure of the superfluid response to 
classical rotation of a boundary as the rotation speed goes 
to zero, which does not map to the quantized rotation of a 
molecule. While this may not be very significant for anal- 
ysis of the molecule-induced local non-superfluid density, 
given that this appears to be independent of the effects 
of molecular rotation^ it should be borne in mind when 
interpreting the meaning of the local superfluid density. 

In order to estimate the local superfluid fraction 
around an impurity molecule, we compute here the den- 
sity distribution of helium atoms participating in per- 
mutation exchanges which result in paths that are long 



compared to the system size. In our previous studies of 
SF 6 4 H eA r and OCS 4 H eA r clusters, 13 i 17 i 26 we found that a 
robust local superfluid distribution in the first shell region 
could be computed for N > 50 by counting the exchange 
paths involving more than six helium atoms. The result 
is robust to the precise value of this cut-off, i.e., similar 
results are obtained with five- or seven-atom exchange 
paths. As pointed out in Ref. UtI this local estimator 
based purely on the exchange path length does not re- 
flect the tensorial nature of the superfluid response that is 
predicted from the linear response definition, Eq. (jSJ . An 
alternative local estimator that is derived from decompo- 
sition of Eq. (J5J into contributions from small cells has 
recently been proposed in Ref. |3?1 This alternative esti- 
mator does possess the correct tensorial nature of the lin- 
ear superfluid response, but unfortunately it is subject to 
extremely large fluctuations. To obtain meaningful val- 
ues within small regions such as the first solvation shell 
of OCS would require an extremely high computational 
effort. (The application of this estimator in Ref. y(J em- 
ployed cylindrical averaging over the length of a linear 
chain of three HCN molecules, which reduces the fluctu- 
ations.) Thus, in order to make a microscopic analysis 
of the effect of the surrounding helium on the moment 
of inertia of the embedded H 2 -OCS complex, we employ 
here the original local estimator of superfluid fraction 
employing long exchange path length. 



III. STRUCTURE AND ROTATIONAL 
CONSTANTS OF FREE OCS-H 2 

We first performed a PIMC calculation on the free H2- 
OCS complex, i.e., without any helium atom present. 
Figure [21 shows the density distribution of the single hy- 
drogen molecule around OCS in the absence of helium. 
The hydrogen molecule is clearly located at and around 
the global minimum of H2-OCS potential. Its average 
position is r = 3.63 ± 0.33 A, 6 = 106.5° ± 8.5°. The 
u ncertainties are evaluated as the fluctuations Aa; = 

\J (x 2 ) — (x) 2 which provide a measure of the widths of 
the single peak in the density distribution. These PIMC 
values of r and 6 extracted from the finite temperature 
density distributions are very similar to the correspond- 
ing values obtained from ground state expectation values, 
namely r = 3.704 and 8 = 105.7ii& The average hydrogen 
position defines an average center of mass of the H2-OCS 
complex. Transforming to the complex body-fixed frame, 
yields the average coordinates r' = 3.51 ± 0.33 A and 
6' = 106.5° ± 8.5°, where r' is the distance from the 
complex center of mass and 0' the polar angle from the 
complex body-fixed z-axis which is defined to be parallel 
to the OCS molecular axis. We see that there is only 
a very slight displacement of the center of mass away 
from the OCS center of mass (~ 0.117 A) as a result of 
complexing with the light hydrogen molecule. 

We compute the moment of inertia tensor of this com- 
plex assuming that the hydrogen is rigidly attached to the 



5 




z (A) 



FIG. 2: Density of H2 in the free OCS-H2 complex, shown in 
the cylindrical coordinates Z and R with origin at the OCS 
center of mass (see Figure 0. 

OCS molecular rotation. Thus, at each imaginary time 
configuration we evaluate the instantaneous moment of 
inertia tensor in the (instantaneous) body-fixed frame of 
the complex which is centered on the center of mass of 
the complex with z-axis parallel to the molecular axis, 
x-axis perpendicular to this and in the OCS-H2 plane, 
and y-axis mutually perpendicular to both x- and z-axes. 
Averaging over the path integral yields the tensor com- 
ponents 

Iij = [Jo]ij + niH 2 {{r 2 $ij - XiXj)), (6) 

where Iq is the gas-phase moment of inertia tensor of 
free OCS, and the coordinates the components of 

the hydrogen position vector r in the body-fixed frame. 
Note that the y-axis is a principal axis of the moment 
of inertia tensor and corresponds to axis c below. We 
then diagonalize this tensor to obtain the three principal 
values for the moments of inertia, namely I a , lb, and 
I c , and their corresponding principal axes a b, and c. 
We obtain the rotational constants of the complex from 
the inverse of the principal moments of inertia, e. g., 
A = h 2 /2I a . 

The resulting moments of inertia, I a , lb, and I c , and 
rotational constants A, B, and C are listed in Table [I] 
(PIMC(l)). The corresponding principal axes are: a in 
the z;r-plane rotated an angle a away from the body-fixed 
z-axis, b perpendicular to this and within the same plane, 
and c perpendicular to the zx-plane. The angle a is also 
listed in Table [I] Figure |3] shows the relation between 
the body fixed and principal axes of the complex. To 
assess the effect of the displacement of the center of mass 
away from the OCS molecule, we have also calculated the 
moment of inertia tensor and corresponding values of A, 
B, C, and a in a body-fixed frame that is centered at 
the OCS center of mass. These values are shown in Ta- 
ble[I]as PIMC(2). The displacement of the center of mass 
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FIG. 3: Schematic of the geometry and principal axes of the 
free OCS-H2 complex, x and z denote the coordinate system 
in the body fixed frame centered on the OCS molecule center 
of mass. The principal axes (a, b) are shown here in this frame. 
The third coordinate y and principal axis c are perpendicular 
to the plane of the paper, with y directed out of the page and 
c into the page. The asterisk denotes the complex center of 
mass, located at 0.117 A away from the molecule center of 
mass on the OCS-H2 axis. Distances are shown to scale. 



away from the OCS axis is seen to have only small effects 
on the rotational constants and the angle a. The PIMC 
results are compared in Table [I] with i) the values ob- 
tained from exact calculations employing close-coupling 
and descendant-weighted diffusion Monte Carlo calcula- 
tions with the same H2-OCS potential^, and ii) recent 
experimental measurements. 

We discuss first the comparison with the exact calcu- 
lations. This comparison allows us to assess the accu- 
racy of the rigid coupling assumption for H2 . The PIMC 
rotational constants are seen to be in excellent agree- 
ment (~ 1%) with the values derived from converged 
energy levels obtained with the close-coupling BOUND 
program^! Consequently we can conclude that the as- 
sumption that H 2 is rigidly coupled to the OCS molec- 
ular rotation is indeed very accurate for OCS-H2. This 
is not the case for the analogous OCS-He complex, as is 
investigated with other methods and discussed in more 
detail in Ref. ^| 
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TABLE I: Moments of inertia (in amu A 2 ) and rotational 
constants (in cm -1 ) of the free OCS-H2 complex, a is the 
angle between the OCS molecular axis and the a principal 
axis. PIMC(l) is obtained with the body-fixed frame cen- 
tered at the true center of mass of the complex. PIMC(2) is 
obtained with the approximation that the body-fixed frame 
is centered at the OCS center of mass. The experimental 
value of a was estimated from the measured intensity ratios 
in Ref. IT2I a — arctan(fj,b/t-i a )- The exact results in col- 
umn 3 derive from rotational energy levels obtained from the 
BOUND program;— (rotation constants and moments of in- 
ertia) and from a ground state average of the moment of iner- 
tia evaluated with exact densities from diffusion Monte Carlo 
calculations (a)^™ 





PIMC(l) 


PIMC(2) 


EXACTi&a 


Expti^ 


la 


21.93(2) 


22.84(2) 


22.29 


22.14 


h 


85.87(2) 


86.39(2) 


84.50 


84.25 


Ic 


108.26(2) 


109.24(2) 


109.28 


109.78 


a 


5.63° 


6.29° 


5.83° 


7.97° 


A 


0.7679(7) 


0.742(3) 


0.7554 


0.7607 


B 


0.1961(1) 


0.1941(5) 


0.1993 


0.1999 


C 


0.1556(1) 


0.1537(5) 


0.1541 


0.1534 



The experimental measurements in the right hand col- 
umn of Table [0 derive from a recent experimental char- 
acterization of the free OCS-H2 complex by infra-red 
spectroscopy^ From spectral fits made to observed tran- 
sitions between low-lying energy levels of the complex, 
it was determined that OCS-H2 is an asymmetric rotor 
at low energies. Table [I] shows that the experimentally 
measured rotational constants (column 3) are in excellent 
agreement (within ~ 0.5%) with the values obtained from 
the exact BOUND calculations (column 2). This con- 
firms that the OCS-H2 interaction potential used herein 
is very accurate, at least in the region around its global 
minima where the single hydrogen is predominantly lo- 
cated. 

We have now confirmed both the accuracy of both the 
OCS-H 2 potential and the rigid coupling approximation 
for the OCS-H 2 complex. Given these two calibrations, 
we can proceed to directly compare the PIMC rotational 
constants (column 1, PIMC(l)) with the experimental 
values (column 3). The agreement here is also excellent 
(~ 0.5 — 3%). The rotational constants A and C, the an- 
gle a, and the free OCS rotational constant (b = 0.20285 
cm -i\3£ can k e used with the ball-and-stick model of 
Hayman et al^ to obtain an independent estimate of 
the structural parameters r' = rt, s and 0' = 8b s - We note 
that the experimentally derived estimates for the ball- 
and-stick parameters, r bs = 3.719 A, 9 bs — 110.8°^ lie 
within the width of the hydrogen density peak in Fig- 
ure El 

The consistently close agreement of the PIMC results 
for all spectroscopic and structural quantities with the 
corresponding experimental quantities therefore allow us 
to conclude that a rigid coupling analysis made with path 
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FIG. 4: Density of H 2 and He in the (OCS-H 2 ) 4 He 6 3 cluster, 
shown in the cylindrical coordinates Z and R with origin at 
the OCS center of mass (see Figure 0. 

integral densities for H2 provides a highly accurate de- 
scription of the free OCS-H2 complex. We note also that 
the quantum averages are consistent with a simple ball- 
and-stick model, although the large size of the quantum 
fluctuations render this model of limited applicability be- 
yond confirming the average structures. 

IV. STRUCTURE AND ROTATIONAL 
CONSTANTS OF OCS-H 2 EMBEDDED IN 
HELIUM 

Figure 0]i shows the density distribution of the sin- 
gle H2 when the OCS-H2 complex is now embedded in 
a cluster of N — 63 helium atoms. The corresponding 
helium density is shown in Figure 2t>. Note that the ori- 
gin in all figures is set to be at the OCS center of mass. 
The hydrogen distribution is seen to change very little 
as a result of the surrounding helium atoms (compare 
Figure 0^ with Figure [2J , except that the single den- 
sity peak is somewhat sharper in the presence of helium. 



7 



This reflects that the presence of the neighboring helium 
atoms makes the binding of the hydrogen to the OCS 
more compact. The average location of the H2 molecule 
is now r = 3.59±0.31 A, 9 = 105.8°±7.5°, very similar to 
its location in the free complex. This location compares 
well with the values r = 3.9 ± 0.2 A, 6 = 111° ± 6° ex- 
tracted from the spectroscopic measurements in Ref. ITol 
The close similarity of the structure of the OCS-H2 com- 
plex with and without solvating helium suggests that the 
rigid coupling of the hydrogen to OCS remains a valid 
approximation when the complex is embedded in helium. 

Looking closely now at the helium distribution, it is 
evident that the first shell structure of this looks very 
similar to the corresponding first shell structure for the 
OCS 4 He64 cluster (see Figure 3 of Ref. 26), with one 
significant difference. This is that the most intense den- 
sity peak at the global minima consists now of only five 
helium atoms, instead of the six helium atoms seen for 
OCS 4 He 6 4. This difference is due to the fact that one 
helium atom is now replaced by a more strongly-bound 
H 2 molecule. Related calculations with larger numbers 
of hydrogen molecules have shown that the six helium 
atoms located at the global minimum in OCS 4 Hejv are 
completely replaced by five hydrogen molecules^ (Sim- 
ilar conclusions have been recently reached from an in- 
tensity analysis of spectral lines in Ref. l35t l 

In order to estimate the rotational constants of the 
OCS-H2 complex inside the larger helium droplets, we 
use the microscopic two-fluid theory of the response of the 
solvating helium density to the molecular rotation^iii 
First, we investigate the local perturbation of helium su- 
perfluidity around the OCS-H2 complex, using the lo- 
cal superfluid estimator based on the exchange length 
of the Feynman paths. As in our previous study for 
the SF 6 - and OCS-doped helium droplets^ we choose 
an exchange length of six or greater as the criterion for 
paths contributing to local superfluidity at any given po- 
sition. The results are relatively insensitive to this cutoff 
length, since the paths are generally split between very 
small paths (2-5 exchanges) and very long paths (~ N 
exchanges). 

The moment of inertia tensor of the complex inside 
the helium droplet is then evaluated by summing the 
rigidly-coupled single hydrogen and the solvating helium 
contributions to the gas-phase moment of inertia tensor 
of OCS: 

hj = [Io]ij + [lH 2 ]ij + [lHe]ij, (7) 

where the hydrogen contribution i# 2 is given by the sec- 
ond term in Eq. |J?J|. We employ here the body-fixed 
frame of the free OCS-H2 complex defined with origin on 
the OCS center of mass (see Section IrTTjl . This is consis- 
tent with the reference frame employed in the analysis of 
the experimental measurements in Ref. llOl 

Within the local two-fluid approach the helium contri- 
bution can be estimated by 




where p ns (r) is the local helium non-superfluid density. 
This is defined as 

Pns(r) = p(r) - p s (f), (9) 

where p{r) is the total number density of helium and 
p s (r) the local superfluid density estimate obtained by 
using the long exchange length criterion. In our previous 
studies of the rotational constants of impurity molecules 
such as OCS inside pure helium-4 clusters, the molecule- 
induced non-superfluid fraction is localized essentially in 
the first shell region^2*i£ Any non-superfluid density in 
the second shell was seen to result from cluster finite 
size effects and decreases as the cluster size increases m& 
However, in this case of the OCS-H2 complex, the pres- 
ence of the hydrogen molecule may induce an additional 
non-superfluid fraction in the second shell region around 
the OCS. We can determine the extent of this effect by 
comparing the local non-superfluid helium density in the 
vicinity of the H2 molecule for (OCS-H2) 4 He63 with the 
local non-superfluid density in the corresponding regions 
for OCS 4 Hc 6 4. Figures |3] (a) and (b) show the non- 
superfluid helium density distributions of A = 63 he- 
lium atoms around the OCS-H2 complex and of N = 64 
helium atoms around the OCS molecule, respectively. 
There is a clear increase in the non-superfluid density 
in the second helium solvation shell around the OCS- 
H2 complex in the vicinity of the H2 molecule. This 
small region of second shell non-superfluid density is ev- 
idently due to the presence of the distinguishable hydro- 
gen molecule nearby. 

In our two-fluid estimate of the effective moment of in- 
ertia, both this additional non-superfluid fraction in the 
second shell and the non-superfluid fraction of the helium 
in the first solvation shell are assumed to be rigidly cou- 
pled to the rotation of the complex. For computational 
convenience, we further assume that the additional non- 
superfluid fraction in the second shell lies in the instan- 
taneous plane including the linear OCS molecule and the 
single hydrogen molecule. This allows the second shell 
contribution to be estimated from integration over the 
density difference between Figures (a) and (b) . Note 
that this provides an azimuthally averaged value for the 
second shell non-superfluid density. More realistically, 
we expect this density to be localized in three dimen- 
sions around the instantaneous H2 azimuthal position. 
Such a three-dimensional distribution will not necessarily 
be accurately described by an azimuthal average. How- 
ever, sampling the relative azimuthal positions of the H2 
and local non-superfluid helium density is computation- 
ally extremely expensive, so we have restricted ourselves 
to the azimuthally averaged approximation here, bearing 
in mind that this may lead to some inaccuracy in the 
second shell contribution to the moment of inertia tensor 
because of the inaccuracy of the underlying mass distri- 
bution. The effective moment of inertia tensor of Eq. Q 
is then diagonalized to yield the principal axes of the hy- 
drogen complex in helium, d, b, and c, and the principal 
moments of inertia tensor / Q , J5, and I c . 
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FIG. 5: Non-superfluid helium density around (a) the OCS- 
H2 complex in (OCS-H2) 4 He63, and around (b) OCS in 
OCS 4 He 64 . The densities are shown in mesh plots with con- 
tours below, in the cylindrical coordinates Z, R defined in 
Figure 

Possible contributions from the superfluid are lim- 
ited by angular momentum and adiabatic following con- 
straints, as discussed in detail in Ref. 0, and are con- 
sequently very small. A recent angular-momcntum- 
consistent calculation for bare OCS in helium has con- 
firmed that the superfluid contribution is negligible in 
this case^& and this is expected to be the same for the 
OCS-H 2 complex. 

The resulting moments of inertia and principal axis 
angle a are listed in Table ITU The average of the three 
PIMC principal moments of inertia is in excellent agree- 
ment with the experimentally derived average (to 4%), 
lying within the statistical error of the path integral es- 
timate. However we note that the PIMC estimates for 
three principal values of the moment of inertia differ from 
their corresponding experimental values I a , h, and I c , 



TABLE II: PIMC moments of inertia (in amu A 2 ), and prin- 
cipal axis angle a of the OCS-H2 complex inside a 4 He63 
cluster, compared with the corresponding experimental val- 
ues obtained from spectroscopic measurements in large helium 
droplets (N ~ 10 3 )f±& 





PIMC 


Expt. 




246(12) 


199(1) 


h 


328(14) 


297(3) 


Ic 


352(14) 


399(6) 


(la +h+ /c)/3 


309(13) 


297(4) 


a 


7.5° 


39(2)° 



by +23%, +10% and —11%, respectively. There are two 
possible contributions to these individual discrepancies. 
First, as discussed above the local superfluid estimator 
based on exchange path length does not reflect the ten- 
sorial nature of the superfluid response, and consequently 
it may not distinguish effectively among the three prin- 
cipal axes a, b and c. Second, as noted above, we have 
approximated the contribution from the second helium 
solvation shell as deriving from a non-superfluid density 
restricted to the OCS-H2 plane. Small inaccuracies in 
the second shell non-superfluid mass distribution result- 
ing from this approximation may give rise to errors in 
the estimated contribution of this shell to the moment 
of inertia tensor. The fact that the individual principal 
moment of inertia values are not achieved at the same 
accuracy as the average value, indicates that the mass 
distribution in the second shell is not predicted with suf- 
ficient accuracy from our simple estimate based on the 
azimuthally averaged non-superfluid density. The differ- 
ence between the experimental value of the principal axis 
angle a and the PIMC estimate of this seen in Table ILTI 
can also be understood to result from this effect, since 
like the individual moments of inertia, a derives from di- 
agonalization of the moment of inertia tensor. Finally, 
any inaccuracy in the individual values of I a , lb, and 
I c will propagate to the estimate of the inertial defect 
A = I c — I a — We obtain a value A = —222 amu 
A 2 . This is negative, in agreement with the experimen- 
tal finding in Ref. 10 and implying a non-planar effective 
complex in helium^i but is larger in magnitude than the 
experimental value A ~ —100 amu A 2 . This discrepancy 
is also consistent with the 10 - 23% inaccuracy in the 
individual moment of inertia values deriving from the ef- 
fect of using an approximate mass distribution from the 
second helium solvation shell. 

In order to assess the importance of the second-shell 
helium contribution to the moment of inertia of the com- 
plex, we show in Table HTlJ the various contributions to 
the moments of inertia, namely, from the H2 density and 
from the non-superfluid helium density in the first and 
second solvation shells, respectively. As can be seen, 
the second-shell helium contribution is quite significant, 
amounting to ~ 18% of the first-shell helium contribu- 
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TABLE III: The moments of inertia (in amu A 2 ) contributions 
from the rigidly-coupled hydrogen (OCS-H2), the first-shell 
non-superfluid of helium (He I), and the second-shell non- 
superfluid of helium (He II) , respectively. 





OCS-H2 


He I 


He II 


Ia 


22.82 


182 


41 


h 


86.10 


219 


23 


Ic 


108.92 


200 


43 


{I a + h + /c)/3 


72.61 


200 


36 



TABLE IV: PIMC rotational constants (in cm -1 ) of the OCS- 
H2 complex inside a 4 Hee3 cluster, compared with the cor- 
responding experimental values obtained from spectroscopic 
measurements in large helium droplets (N ~ 10 3 )^- 





PIMC 


Expt. 


A 


0.069(4) 


0.0847(4) 


B 


0.051(3) 


0.0567(6) 


C 


0.048(3) 


0.04226(6) 


(A + B + C)/3 


0.056(3) 


0.04948(6) 



tion to the average moment of inertia of the complex. 
We find that if the contributions from the second solva- 
tion shell are omitted (OCS-H2 column plus He I column 
in Table ITTTfl . the average moment of inertia differs from 
the experimental quantity by a greater amount, namely 
—8%, while the individual principal moments of inertia 
differ now from experiment by +3%, +2% and —23%, re- 
spectively. The fact that the average moment of inertia 
is now significantly lower and in poorer agreement with 
the corresponding experimental value confirms the need 
to take the second shell non-superfluid contribution into 
account. 

Since we calculate directly the moment of inertia ten- 
sor of the OCS-H2 complex, a comparison of the prin- 
cipal moments of inertia with corresponding experimen- 
tal values is the best way for us to calibrate our accu- 
racy. Assessment of the resulting rotational constants 
is also possible, but is complicated by the need to in- 
vert the principal values and the consequent increase in 
error. Table IIVI shows the comparison of the PIMC ro- 
tational constants and the corresponding experimental 
values. Even though the average rotational constant is 
less accurate than the average principal moment of iner- 
tia, the PIMC-based estimate for the average rotational 
constant is still in good agreement with the experimental 
value, to ~ 15%. Given the approximate description of 
the second shell helium contribution, this would appear 
to be excellent agreement. 

V. DISCUSSION 

We have presented microscopic quantum calculations 
of the rotational constants of OCS-H2 both as a free com- 
plex and when embedded in helium droplets. The analy- 



sis was based on path integral calculations, and employed 
an assumption of rigid coupling of the H2 motion to the 
OCS rotation, as well as the microscopic local two-fluid 
theory to evaluate the helium contribution to the mo- 
ment of inertia. The OCS-H2 complex is found to be an 
asymmetric rotor, and to possess a very similar structure 
when free and when embedded in helium. The complex 
is seen to induce a local non-superfluid helium density in 
the second solvation shell, in the immediate vicinity of 
the H2 molecule, as well as in the first shell. The overall 
effect of the solvating helium is to change the complex 
geometry as measured spectroscopically, from planar to 
non-planar, reflecting an additional moment of inertia 
contribution from the solvating helium. The calculated 
moments of inertia derived from the two-fluid estimates 
and corresponding rotational constants are in excellent 
agreement with experimentally measured values for the 
free OCS-H2 complex, with individual values within 3%. 
For the OCS-H2 complex embedded in a helium cluster, 
the average moment of inertia is in excellent agreement 
with the corresponding experimental average, within 4%, 
when the contribution from the non-superfluid density in- 
duced by the H2 molecule in the second solvation shell 
of helium is taken into account. The individual values of 
the moments of inertia have larger deviations from exper- 
iment, 10 — 23%, but are still in good agreement. These 
individual deviations, as also the deviation in angle a 
and in magnitude of inertial defect A may derive from 
inaccuracies in the components of the moment of iner- 
tia tensor resulting from treating the second shell helium 
contribution as azimuthally symmetric. 

The two-fluid analysis made here provides a physically 
consistent explanation of the anomalous effective masses 
assigned to the H9 molecule in the extended shell model 
proposed in Ref. lloL In that model, an extension of 
the earlier donut model for OCS in helium droplets, 38 
the increased moment of inertia was assigned to con- 
tributions from H2 and from non-superfluid helium in 
the first solvation shell only. The latter was estimated 
from a close packing model of our PIMC helium density 
distributions ^£ with effective helium masses of 0.55 amu 
assigned to approximate the theoretical non-superfluid 
density. The contribution from the H2 molecule was 
described with an adjustable effective mass parameter, 
whose best fit value to experiment was 10 amuM^ This 
large increase in effective mass of H2 was attributed in 
Ref. ^3 t° the substitution of one indistinguishable he- 
lium atom by a distinguishable H2 molecule, without any 
associated change in the surrounding helium density dis- 
tribution. The analysis made here provides a simpler and 
more physical explanation for the anomalously large ef- 
fective mass of H2, namely, that this necessarily results 
when the HVinduced second shell non-superfluid density 
is neglected and the moment of inertia increment is as- 
sumed to result only from the first solvation shell density 
components. The fitted value of the H2 effective mass 
is larger than its natural mass because this adjustable 
parameter must absorb the second shell contributions in 
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the extended shell model of Ref. It is therefore iden- 
tified as an artifact of the model that would be removed 
if the model were extended to include the second shell 
helium contribution. 

A rigid coupling of the hydrogen molecule to rotation 
of a heavy molecule such as OCS seen here is not unex- 
pected. Related studies have shown that this description 
becomes less accurate for higher total angular momentum 
states of the complex^ It is useful to compare this be- 
havior of the OCS-H 2 complex with that of its cousin, 
the OCS-He complex, which is considerably less rigid 
(see, e.g., the comparative analysis in Ref. Ilq) . Despite 
having similar topology of interaction potentials, the H2 
molecule has considerably stronger binding to OCS than 
does He. This results in a more sharply modulated struc- 
ture, less derealization, and more rigidity in coupling of 
its density to the OCS motion, reflecting the presence of 
low lying excitations of a rigid complex. Nevertheless, 
it is important to realize that the hydrogen component 
is still a highly quantum element and this becomes very 
evident for larger M (see below). 

The excellent agreement with a rigid coupling model 
for the free complex suggests that similar rigid cou- 
pling, or more general quasi-adiabatic coupling; 16 ; 39 de- 
scriptions might provide good approximations for the 
analysis of OCS(H2)m, both free and embedded in he- 
lium droplets. However, at larger sizes, in particular for 
M > 11, we have recently shown that the hydrogen com- 
ponent possesses an anisotropic superfluid stated This 
state is by definition characterized by a lack of rigid re- 
sponse to rotation about the molecular axis. The first 
solvation shell is complete at M = 17, and at interme- 
diate M values one can identify precursors to the super- 
fluid state that possess significant permutation exchanges 
within certain segments of the local solvation structured 
Thus for intermediate sizes one expects the values of ro- 
tational constants to reflect a balance between rigidity 
due to the strong OCS-H2 and H2-H2 interactions, and 
lack of rigidity due to onset of these permutation ex- 



changes. The intermediate-size hydrogen complexes with 
OCS provide a rich opportunity to analyze the detailed 
evolution from full rigidity to anisotropic rigidity and su- 
perfluidity. 

The high accuracy of this two- fluid analysis for the mo- 
ments of inertia and rotational constants of the OCS-H2 
complex when embedded in helium provides confirmation 
of the accuracy of the underlying local two-fluid theory. 
This analysis of moments of inertia and rotational con- 
stants for OCS-H2 inside a helium droplet provides the 
first microscopic quantum calculation of rotational con- 
stants for a van der Waals complex inside helium. The 
average values are found to be highly accurate, and the 
slightly larger deviations for the individual values un- 
derstandable in terms of simplifications of the helium 
non-superfluid density made here to allow computational 
tractability. The present analysis thereby extends the 
previous findings of highly accurate two-fluid calculation 
of rotational constants for the isolated molecules SFg and 
OCS in helium cluster a 13 ' 17 to embedded molecular clus- 
ters. In the case of these heavy molecules and clusters for 
which high quality interaction potentials are available to 
enable meaningful microscopic calculations to be made, 
the two-fluid theory thus appears to provide an accurate 
microscopic description of response of a solvating super- 
fluid helium droplet to the molecular rotation. 
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